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ABSTRACT 


A 


Nixing  lengths  ere  formulated  for  turbulent  boundary 
layers  over  arbitrarily  rough  surfaces  from  both  the 
van  Driest  and  fro*  the  Rotta  procedures.  The  associated 
eddy  viscosities  and  the  associated  turbulent  kinetic 
energies  and  turbulent  dissipation  rates  are  also  formulated. 
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ADMINISTRATIVE  INFORMATION 

The  work  In  this  report  was  performed  at  the  David  U.  Taylor  Naval  Ship  Research 
and  Development  Center.  It  ms  authorised  and  funded  by  the  General  Hydromechanics 
Program  of  the  Naval  Sea  Systems  Command  under  Program  Element  M135N,  Task  Area 
SR0230101.  and  Work  Unit  1542-070. 


INTRODUCTION 

The  development  of  turbulent  boundary  layers  over  arbitrarily  rough  surfaces 

nay  be  predicted  by  mixing-length  formulations  close  to  the  boundary  wall.  Existing 

calculation  methods  for  smooth  surfaces  using  the  Prendtl-ven  Driest  well  nixing 

length  nay  be  adapted  to  arbitrarily  rough  eurfaces.  Such  methods  for  smooth  eur- 

faces  are  the  method  of  Patankar  and  Spalding.  the  method  of  Cebecl  and  Smith,'  or 

j  u 

Its  subsequent  modification  by  Nituch  at  al.  and  the  method  of  Huang  et  al . 

Wall  mixing  lengths,  or  rather  the  aseoclated  eddy  viscosities,  sre  also 
required  for  heat  or  mass  transfer  methods  using  turbulent  Prsndtl  numbers  ot  turbu¬ 
lent  Schmidt  numbers  to  predict  the  turbulent  dlf fuslvltlee  of  heat  or  mass. 
Furthermore,  mixing  lengths  or  eddy  viscosities  near  the  wall  auiy  he  required  as 
boundary  conditions  to  the  k*f  transport  equations  for  improved  predictions  ol  the 
development  of  turbulent  boundary  layers.6  tn  general,  the  mixing  length  increases 
with  distance  away  from  the  well  until  e  limiting  value  Is  reached  which  is  compat¬ 
ible  to  the  mixing  length  for  the  outer  region. 

Initially,  In  1950,  Rotta7  considered  the  Prendtl  wall  mixing  length  for 
turbulent  shear  stress  to  epply  only  outside  the  laminar  sublayer  and  was  able  to 
correlate  the  laminar  sublayer  thickness  with  the  roughness  drag  character  teat  ion 
function.  However,  subsequent  turbulence  measurements  showed  that  turbulent  shear 
stresses  start  at  the  wall.  The  concept  of  e  purely  laminar  sublayer  was  discarded 
and  the  laminar  sublayer  was  renamed  the  viscous  sublayer.  For  sswoth  surfaces,  the 


*A  complete  listing  of  references  Is  given  on  page  31. 
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van  Driest  Modification  to  the  Prandtl  wall  mixing- length  formulation  provided  a 

a 

turbulent  shear  stress  which  started  at  the  wall.  Rotta  then  adapted  the 

van  Dr  lest ^Prandtl  nixing  length  for  moo eh  surfaces  to  rough  surfaces  by  shifting 

the  position  of  the  reference  wall  an  appropriate  distance.  Rotta  presented  tie 

masericai  results  in  graphical  form  for  two  well-known  roughnesses:  the  Nlkuradse 

sand-grain  roughness  and  the  Col ebrook- White  engineering  roughness,  hater  Cebet i 

and  Snlth^  provided  en  analytical  fit  to  the  wall  shift  for  the  sand-grain  roughness 

9 

which  wee  then  used  by  Cebeci  end  Chang  for  calculations  of  turbulent  boundary 
layers  over  rough  surfaces. 

Tor  arbitrarily  rough  surfaces,  the  van  Driest  factor  is  correlated  in  this 
paper  with  the  roughness  drag  characterisation  function  until  a  limiting  value  of 
sero  is  reached.  This  Is  considered  to  represent  the  beginning  ol  the  fully  rough 
regime.  Then,  for  the  fully  rough  regime,  the  nixing  length  is  considered  to  assume 
an  initial  value  at  the  wall  in  accordance  with  Rotta'*  original  (I9J0)  analysis. 

For  purposes  of  comparison,  the  second  Rotta  analysis  of  a  wall  shift  is 
generalised  to  produce  numerical  result*  for  the  drag  eharartet last  Ion  of  arbitrary 
roughness,  tn  Appendix  A,  an  explicit  conversion  of  n lx li*- length*  to  eddy  viscos¬ 
ities  is  derived.  In  Appendix  I,  the  sssoctated  turbulent  kinetic  energies  (f. )  and 
turbulent  dissipation  rate  (<)  are  developed. 

CHARACTERIZATION  OF  ROtXHXfcSh  DRAT. 

A  brief  review  is  presented  of  pertinent  feature*  the  velocity  similarity 
law*  and  the  associated  drag  character  i  tat  ton*  c(  rough  surface*. 

Close  to  the  wall,  the  logariltmic  velocity  law  (or  an  arbitrarily  rough 

surface,  defined  by  a  sufficient  number  of  length  factors  k,  ij,  i, . i«  given 

tn  what  nay  be  catted  the  Re  mo  ids- number  mode  a* 

—  «  A  in  y*  ♦  t.  |k*,TJ  <1  * 

Mt  1 

and.  In  what  nay  be  called  the  relat tee- roughness  mode,  as 

--  *  A  In  l  *  «  |k*,T|  (2) 

u  ^  k  r 

where 


2 


(3) 


•  Ej  ♦  A  In  k* 

Mere  u  •  »tr««mlw  velocity  component 
•  sheer  velocity,  u?  • 
i  •  well  cheer  stress 
w  *  density  of  fluid 
* 

y  *  u  v/v 

T 

y  •  no reel  distance  f ran  well 
v  •  kinematic  viscosity  of  fluid 

t  •  texture  of  roughness  conf igurat ion,  T  *  k/ij,  kj/kj  ... 

•  s 

k  •  roughness  Reynolds  nunher,  k  •  «(k/v 

Cither  1.  or  I  ney  he  considered  to  be  a  roughness  dreg  character irai ion 
*  *  • 
function.  Roth  Bj  end  B,  ate  functions  of  roughness  Reynold#  nunher  k  and 

texture  7. 

The  usual  roughness  regimes  ere 
1.  Mydveul Icaily  snoot h 

R,  *  R,  •  constant 
1  1  ,s 

t  *1.  •  A  to  k 

r  I .» 

*.  Intermediate  rowghress 

Roth  R|  and  lf  vary  with  k  for  the  sane  1 

).  Fully  rough 


•t  *  lf  | T|  •  constant 

,  <*» 

•j  *  Rf  |  T|  -  A  In  k 

Another  dreg  character  I ration  Is  given  by  AR  which  represents  a  deviation  from 


3 


smooth  condition*  or 


AB  |k‘,T!  -  B.  I k* , T J  -  B. 

1  1 ,  s 


(*>) 


As  defined,  AB  Is  always  negative  for  rough  surfaces.  Sikuradse  used  the  B  - 

11  r 
characterization  while  Hama  preferred  the  AB  (actually  -AB)  characterization. 

fo  experimentally  determine  a  characterization  tor  a  specific  arbitrarily  rough 
surface,  there  are  various  procedures  available.  The  direct  procedure  requires 
velocity  measurements  close  to  the  wall  to  define  a  logarithmic  law  as  well  as  the 
measurement  of  the  wall  shear  stress.  Simpler,  Indirect  procedures  may  he  used 

Involving  measurement  of  the  average  velocity  of  pipe  flow,*0  the  total  drag  of  a 

1  -  13 

flat  plate,  or  the  torque  of  a  rotating  disk. 

MIXING- LENGTH  THEORY 

According  to  the  Brandt  1  m lx ing- length  theory,  the  turbulent  shear  stress  ( '  ) 
to  related  to  the  velocity  gradient  by 


t  .  2  i  du 


(7) 


where  I-  U  the  mixing  length  which.  In  general,  is  not  constant.  but  >•  'unction 
.to  po«  it  ion  in  the  flow  field. 

The  Cota!  shear  stress  '  at  a  point  In  the  flow  :  »«Td  i  •.  the  mi-,  of  laminar 
*pd  turbulent  -  shear  drr««r<  or 

t 


du  , 2  /  du 

*  *  I 


,  the  Newtonian  Law  of  Viscosity. 

For  a  shear  layer  with  zero  longitudinal  pressure  gradient,  '  -  ~  Close  to 
wall.  With  this  approximation,  Equation  (8)  mav  be  written  in  a  nond i -a*ns iona 1  torr 


(<?) 


4 


* 

where  t  is  a  nondi»en« tonal  aUln^  length,  l  *  u  t/vf  and  u  it  a  nondlaens tonal 
*  T 
veloc ity ,  u  •  u/u? * 

a  a 

Solving  for  du  /dy  aa  a  quadratic  expression  and  integrating  from  the  wall, 
a  a 

where  y  »  0  and  u  •  0,  result «  In  a  velocity  profile 


a 

u 


dy 


l  ♦  /t  +  (2iV 


(10) 


Close  to  the  wall  Prandtl  proposed  a  linear  variation  of  mixing  length  with 
distance  fro*  the  wall 


» y 


(11) 


where  »*  l#  the  von  Kansan  constant.  Use  of  this  relation.  In  Equation  (10),  leads 

to  a  logarithmic  veloc itv  law  with  A  *  !/• .  However,  an  erroneous  value  of  Bj 

results.  Furthermore,  In  so  far  as  the  logarithmic  velocity  law  does  not  hold 

right  up  to  the  wall,  so  the  Prandtl  wall  Mixing  length  does  not  hold.  To  remedy 

#  «  *  .  * 
this,  a  modification  function  M(*.  “JV'y  )  Is  needed  as  a  function  of  y  and  <  where 
*  * 

<  is  an  additional  nondiaens tonal  length,  *  *  u.«/v.  Furthensore,  M  should  equal 

*  • 
zero  at  v  *0  and  equal  unity  at  y  •* 

14 

For  smooth  surfaces,  van  Driest  proposed  the  representation 


M 


(12) 


where  A  Is  the  A  associated  with  van  Driest, 
v 

Other  modification  functions  stay  he  formulated  as,  for  example,  the  hyperbolic 
tangent  function  alluded  to  by  Patel ^ 


where  A.  is  the  A  used  here, 
n 


M  - 


tanh 


(13) 


i 


Theo- 

*1 


Why  H  should  equal  taro  at  y  •  0  requires  sone  further  discussion, 
ret  leal  investigations^  indicate  that  the  eddy  viscosity  v  thould  varv  with  y  1  or 

*4  *  t 

y  at  y  «  0.  Mow  Equation  (A. 6)  in  Appendix  A  relates  the  eddy  viscosity  to  thc- 
mixing  length.  A  Haclaurin  expansion  gives 


v 

J 

V 


2 


<1M 


or 


v 


2*2  y*2  M2  ♦  ... 


A  nonzero  value  of  M  at  y  •  0  tuch  a*  Hj  would  then  rcault  in 


2<2  H,2  y*2  ♦  . . . 


(IS) 


which  is  not  acceptable. 

A  Mnclaurtn  expansion  of  the  van  Driest  H  gives 


e 


v 


>nJ 


v 


which  1 9  acceptable. 


6 


UIimIm,  for  th*  hypo  r  bo  lie  t«o|«at  M 


v 

t 

v 


♦ 


(18) 


EXISTING  MIXING- LENCTH  METHODS  FOR  ROUGH  SURFACES 
Som  port  loom  existing  six  lag-length  methods  for  rough  ourfocoo  oro  now 
critically  examined . 

ROTTA  1 

An  extension  of  the  Frendtl  wall  alxlng  length  net hod  to  rough  ourfocoo, 
initially  proposed  by  Rotta, ^  assumed  a  laminar  oublayor  thtcknoao  that  would  de¬ 
er  oaao  with  roughfiooo  or 


l*  -  *(y*-y^  (19) 

*  * 

whore  y^  la  the  laalnar  oublayor  thlcknooo  and  y^  ■  uTyL/v.  The  Introduction  of  i 

In  Equation  (10)  resulted  In  an  Integration  In  oloaentary  functiono  and  a  relation 

a 

between  y^  and  B^ 


y*  -  Bj  ♦  i  (1-  In  4*)  (20) 

* 

With  increaaing  roughness,  the  laalnar  sublayer  finally  vanishes,  yL  •  0,  and 
Bj  •  -  1/*  (l-ln  4k).  This  is  now  considered  to  be  the  beginning  of  the  fully 
rough  regime  and  the  end  of  the  intermediate  roughness  regime. 

For  the  fully  rough  regime,  Rotta  aoouaed  an  Initial  alxlng  length  at  the 

wall 


a  *  a  „ 

t  -  l  +  pey  (21) 

w  1 

The  introduction  of  this  alxlng  length  for  the  fully  rough  regime  into  the 
velocity  profile.  Equation  (10),  results  in  an  analytical  solution. 


(22) 


Ac  higher  values  of  y*,  this  relation  reduces  to  the  usual  log  law,  Equation  (1),  so 
that 


B 


1 


(23) 


* 

This  relates  B^  and  therefore,  the  velocity  profile.  Equation  (10),  Is  a 

* 

function  of  y  and 


ROTTA  II 

8 

Rotta  seems  to  have  abandoned  the  laminar  sublayer  approach  in  analyzing 
rough  surfaces  after  the  whole  concept  of  a  purely  laminar  flow  next  to  the  wall 
became  untenable.  Since  there  are  turbulent  stresses  even  in  the  laminar  sublayer, 
the  sublayer  has  been  renamed  the  viscous  sublayer. 

Rotta  extended  the  van  Driest  formulation  for  smooth  surfaces  to  rough  surfaces 
by  adding  a  length  Ay  to  the  y  coordinate  such  that 


*  *  1 

*(y  +Ay  )  ( 1  -exp 


l 

1  v,s 


(24) 


8 


*  * 

where  A  is  the  velue  of  A  for  saooth  surf sees.  When  this  mixing  length  is  in- 
vi»  v 

corporsted  into  the  velocity  profile,  Equation  (10),  and  related  to  the  logarithmic 

*  * 

law.  Equation  (1),  Ay  becomes  a  function  of  B.(k  ,T].  Instead  of  solving  the 
*  *  1  * 

general  case  of  Ay  ■  flB.J,  Rotts  solved  Ay  ■  f(k  )  for  the  Nikuradse  sand-grain 
*  e  * 

roughness  and  Ay  ■  f | k  ]  for  the  Colebrook-Whlte  engineering  roughness  and  presented 

2 

the  results  graphically.  Cebecl  and  Smith  fitted  the  results  for  the  Nikuradse 

o 

sand-grain  roughness  to  an  empirical  formula.  Subsequently,  Cebeci  and  Chang  used 
this  formula  in  boundary-layer  calculations. 

*  * 

The  Rotta  formulation  leads  to  an  initial  value  of  £  ,  namely,  l  at  the  wall 
*  w 
(for  y  «0)  or 


Hence  •  f[Bjl  -  f [ k  ,T]  for  all  roughness  regimes. 


(25) 


DAHM 

Dahm*6  proposed  a  differential  equation  for  the  mixing  lengths  of  rough  sur¬ 
faces  based  on  empirical  considerations,  which,  for  zero  mass  injection,  becomes 


* 

dl 


*  *  * 

0.4  y  -  (i  -*w) 

TT783 


(26) 


Examination  indicates  that  this  equation  is  a  linear  differential  equation  with  a 
solution  given  by 


t*  -  0.4  y*  +  l*  +  (0.4)  (11.83) 


(— r1 - A 

\ey  /n-83  ) 


(27) 


9 


i 


It  uy  be  noted  that  this  formula  resembles  that  of  fcotta,  Equation  (21),  for  the 

ASS  A 

fully  rough  regime  (x»0.4;  also  at  y  *0,  and  t  “1.^).  In  actual  use.  values  of  t w 
would  have  to  be  correlated  with  drag  characterisation  (B^). 


PROPOSED  NEAR-WALL  MIXINC  LENGTHS  FOR  ROUGH  SURFACES 


GENERAL 


Equating  the  velocity  profiles  obtained  from  mixing  lengths.  Equation  (10),  and 
the  logarlthailc  law.  Equation  (1),  produces 


and  with 


Vi 


0  1  +  /l+(2l*)2 


*  1 

dy  -  -  In  y 


ry  * 

M  ^ 

j .  V 


rl  ry  r  -I 

-  f  - } - dy*  +  f  r -  2  -  --  ~  dy* 

J0  1  +  /+<2iV  J1  [_1  +  ^+(2l  )2  Ky  _ 


B1  -  h  +  X2 


where  y*  is  a  sufficiently  large  value  of  y  so  that  the  second  integrand  of  Equation 


(30),  dl2/dy  ,  becomes  negligible. 


INTERMEDIATE  ROUGHNESS  REGIME 

With  the  Prandtl-van  Driest  mixing  length. 


L-  ex 
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in  Equation  (30),  the  van  Driest  factor  A  becomes  a  function  of  B.  and,  hence,  a 

*  v  *  *  1 
function  of  roughness  factors,  k  and  T  or  A  ■  f ( )  >  f(k  ,T|. 

*  vx 

Numerically  then,  A  from  the  preceding  equation  decreases  with  decreasing 
v  *  *  * 

values  of  B.  until  the  lialt  A  ■  0  Is  reached.  Then,  1  ■  <y  .  This  Is  the  saae 

Halt  as  Rotta  I  where  the  laalnar  sublayer  disappears. 

*  * 

The  range  of  A  from  Its  saooth  value  of  A  to  Its  zero  Halt,  may  be 
v  v » 8  *  * 

considered  to  cover  the  Interned  late  roughness  reglne,  at  Ay  ^  >  A^  >  0. 

FULLY  ROUGH  REGIME 

* 

The  original  Rotta  I  proposal  of  an  initial  nixing  length  1  at  the  wall  for 
the  fully  rough  reglne  will  now  be  considered  appropriate.  Accordingly,  Equations 
(11).  (22),  and  (23)  apply. 

GENERALIZED  ROTTA  II 

In  order  to  make  comparisons,  the  Rotta  II  method  of  an  additional  length  Ay 

Is  now  generalized.  The  Rotta  II  mixing  length  for  rough  surfaces,  Equation  (24), 

* 

is  substituted  Into  either  Equation  (28)  or  (30)  for  B.  so  that  Ay  becones  a 

*  *  1 
function  of  ,  so  that  Ay  «  f ( B^ |  -  f(k  ,Tj. 

DISCUSSION 

Two  methods  have  been  described  for  the  prediction  of  mixing  lengths  for 
arbitrarily  rough  surfaces.  There  Is  the  Rotta  II  method  which  is  generalized  in 
terms  of  the  roughness  characterization  function  and  applies  to  both  the  intermediate 
and  fully  rough  regimes.  The  Rotta  II  method  extends  the  van  Driest  formula  for 
smooth  surfaces  to  rough  surfaces  by  means  of  a  normal  distance  parameter. 

The  other  method  described  consists  of  an  extension  of  the  van  Driest  formula 
to  the  intermediate  roughness  regime;  the  van  Driest  factor  is  correlated  to  the 
roughness  characterization  function.  For  the  fully  rough  regime  the  Rotta  I  method 
is  to  be  used  In  which  a  wall  mixing  length  is  related  to  the  roughness  character¬ 
ization  function.  The  relative  merits  of  each  method  are  now  discussed. 

For  the  intermediate  roughness  regime,  the  Rotta  II  method  leads  to  an  initial 
mixing  length,  Equation  (2S),  at  the  wall  while  the  van  Driest  factor  method  gives 


11 


*  zero  mixing  length.  A  nonzero  nixing  length  at  the  wall  inpliea  the  existence  of 
a  nonzero  turbulent  ahear  atreas  and  a  nonzero  eddy  viacoaity  at  the  wall  aa  la 
deaonetrated  next. 

Now  the  turbulent  ahear  atreaa  near  the  wall.  Equation  (7),  nondinenaionallzed 
to 


*  * 

becomes,  by  use  of  Equation  (10)  for  du  / dy 


r  * 

_ n _ I 

_l+/+<2tVj 


then,  at  the  wall  y  *  0, 


l+/+(2t*)2 


03) 


(34) 


(35) 


Consequently,  at  the  wall,  the  Rotta  II  method  has  a  wall  value  of  turbulent  shear 
stress,  but  the  van  Driest  factor  method  has  none. 

Also  for  the  eddy  viscosity.  Equation  (A. 6)  gives 


\J' 


*  2 

l+(2lw) 


1 

2 


(36) 


Again,  the  Rotta  II  method  has  an  initial  value  of  eddy  viscosity  at  the  wall  while 

the  van  Driest  factor  method  has  none.  In  fact,  the  van  Driest  factor  method 

applied  to  the  intermediate  rough  regime  regains  the  variation  of  eddy  viscosity 
*& 

with  y  at  the  wall. 
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For  Che  Intermediate  roughness  regime,  an  initial  value  of  mixing  length, 
turbulent  shear  stress,  and  eddy  viscosity  at  the  wall  may  all  be  considered 
object lonable. 

NUMERICAL  RESULTS 

In  determlng  B.  from  an  Integration  of  mixing  lengths  by  means  of  Equation 

(30),  It  la  necessary  to  specify  a  limiting  value  of  y  so  that  the  Integrand  of  I. 

*  *  2 

becomes  practically  zero.  As  an  example,  the  integrand  v  (dl./dy  )  Is  plotted  In 

a  2 

Figure  l  for  the  smooth  case  X  *26  and  for  the  case  of  the  beginning  of  the  fully 

#  v 

rough  regime  X  ■  0  with  <  ■  0.4  (the  same  case  as  for  ROTTA  II).  To  compensate  for 
v  * 
the  logarithmic  abscissa,  the  ordinate  is  multiplied  by  y  so  that  the  area  under 

*  *  *  *  * 

the  curve  represents  the  Integral  I,.  Then,  dl,/dy  ■  y  (dl,/dy  )  (d  f.ny  / dy  ). 

Ac  values  of  y  above  about  200  the  integrand  is  close  to  zero. 

* 

In  Figure  2,  the  van  Driest  factor  X^  is  correlated  with  drag  characterization 

function  B.  for  the  intermediate  roughness  regime  by  means  of  Equation  (33).  This 
1  * 
accommodates  any  arbitrary  roughness.  The  limiting  value  Is  Bj  -  >1.323  for  >v  *  0. 

Smaller  values  of  B.  are  then  In  the  fully  rough  regime. 

1  * 

The  wall  mixing  length  t  for  the  fully  rough  regime  is  plotted  against  the 
drag  characterization  function  In  Figure  3  In  accordance  with  the  Rotta  1 
relation.  Equation  (23). 

* 

In  Figures  4  and  5  the  Rotta  factor  Ay  Is  plotted  against  the  drag  character¬ 
ization  function  B^,  as  a  result  of  solving  Equations  (24)  and  (30).  This  corre¬ 
lation  applies  to  any  arbitrary  rough  surface  and,  consequently,  includes  the 

Nlkuradse  sand-grain  roughness  and  the  Colebrook-Whlte  roughness. 

* 

A  comparison  is  shown  in  Figure  6  for  the  mixing  lpngtb  •  at  the  wall  as  a 

function  of  the  roughness  characterization  function  B.  .  It  should  be  noted  that  the 

*  1 

Rotta  II  formulation  has  values  of  i  in  the  intermediate  roughness  regime  in 
accordance  with  Equation  (25).  On  the  other  hand,  the  van  Driest  formulation  has  a 

zero  value.  For  the  fully  rough  regime  the  Rotta  I  values  are  compared  to  the 

* 

Rotta  II  values  of  l  . 

w  *  * 

A  comparison  of  the  variation  of  mixing  lengths  ■  with  normal  distance  y  is 

displayed  in  Figure  7  for  various  formulations  on  the  basis  of  equal  values  of  Bj . 
The  smooth  case  is  shown  for  a  van  Driest  factor  of  26.  As  an  example  of  the 
intermediate  roughness  regime,  a  comparison  is  shown  between  the  van  Driest  and 


Rotta  II  formulations  for  Bj  -  2.  It  is  to  be  noted  that  the  Rotta  II  value  of  the 

nixing  length  la  larger  close  to  the  wall.  This  is  a  result  of  having  an  initial 

* 

value  l  at  the  wall.  Mixing  lengths  at  the  border  between  the  interned  late 
"  * 
roughness  and  fully  rough  regiaes  are  also  compared.  Here  X^  ■  0  and  *  -1.125. 

For  the  fully  rough  reglae.  an  example  comparison  is  shown  for  B^  ■  -6.  It  is  to 

*  * 

be  further  noted  that  for  large  values  of  y  all  the  values  of  i  tend  to  converge 

*  * 

ro  the  original  Prandtl  wall  nixing  length  of  1  •  <y  as  expected. 

SUMMARY 

For  arbitrarily  rough  surfaces,  two  procedures  have  been  developed.  First  is 

the  extension  of  the  Prandtl-van  Driest  formulation.  Equation  (12),  to  the  lnter- 

* 

mediate  roughness  regime.  Here  the  van  Driest  factor  is  correlated  with  drag 

characterization  B.  in  Figure  2.  For  the  fully  rough  regime,  the  Rotta  1  procedure 
1  * 
miy  be  used  as  given  by  Equation  (23)  which  is  plotted  in  Figure  3  as  l against  . 

Second,  there  is  the  Rotta  II  correlation.  Equation  (2A),  developed  for  arbitrarily 

a 

rough  surfaces.  Figure  A  shows  the  correlation  of  Ay  with  B^  for  the  intermediate 
roughness  regime  and  Figure  5  shows  that  for  the  fully  rough  reg<me. 
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Figure*  1  -  Variation  of  Integrand  v  I  — ^  I  with 
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NONOIMENSIONAL  MIXING  LENGTH. 


APPENDIX  A 

CONVERSION  OF  MIXING  LENGTHS  TO  EDDY  VISCOSITIES 


The  concept  of  eddy  viscosity  v  (more  properly  eddy  kinematic  viscosity)  was 
originally  proposed  during  the  last  century  by  Boussinesq  for  turbulent  flows  as  an 
analogy  to  the  usual  viscosity  for  laminar  flows  or 


T  t  du 

P  Vt  dy 


(A.  1 ) 


Equating  this  to  the  turbulent  shear  stress  given  in  terms  of  mixing  lengths. 
Equation  (7),  produces  the  well  known  expression 


v 


t 


l1 

dy 


(A. 2) 


Note  that  eddy  viscosity  unlike  laminar  viscosity  is  not  a  property  of  the  fluid, 
but  depends  upon  its  position  in  the  flow. 

For  various  reasons  it  may  be  more  desirable  to  relate  eddy  viscosity  (v^)  to 
mixing  length  without  the  presence  of  velocity  gradient  du/dy. 

In  general,  shear  stress  (t)  has  laminar  and  turbulent  contributions  such  that 


jr 

P 


(v+vt) 


du 

dy 


(A. 3) 


Close  to  the  wall  T  ~  Tw  and  then,  nondimensionally , 


Also  from  Equation  (10) 


(A. 4) 


* 

du  _  _ 2 _ 

dy  1  +  /+( 2d*)2 


(A. 5) 
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ttoonzio  jus  aum-mr  num 


*  * 

Equating  the  two  expressions  for  du  /dy  results  in 

Vt  1  /  *  2  1 

^  j  /l+(2i l  )  -  ~  (A. 6) 

This  is  the  desired  relation  for  converting  mixing  lengths  to  eddy  viscosities. 

*  *  ★ 

At  large  values  of  l  ,  (v^/v)  -*•  il  ■*  <y  or  -*■  u^S.  -*■  Ku^y. 

For  the  van  Driest  formulation  for  arbitrarily  rough  surfaces.  Figure  A.l  shows 
the  variation  of  eddy  viscosity  ratio  V  /v  with  normal  distance  y  .  The  smooth  case 
is  shown  as  well  as  the  case  for  the  boundary  between  the  intermediate  and  fully 
rough  regimes.  Also,  as  an  example  of  the  fully  rough  regime  v  /v  is  shown  for  B.  = 

it  a  ^  I 

-6.  At  large  values  of  y  ,  (v  /v)  -*  <y  as  seen  in  Figure  A.l. 

c  * 

It  is  to  be  noted  that  the  increase  of  \)  /\>  with  y  presented  here  applies  only 

to  the  boundary-layer  region  next  to  the  wall.  A  limiting  value  is  reached  which 

2 

corresponds  to  outer-region  values  such  as  those  given  by  Cebeci  and  Smith. 


2U 


Figure 


APPENDIX  B 


NEAR-WALL  VALUES  OF  TURBULENT  KINETIC  ENERGY  (fe) 
AND  TURBULENT  DISSIPATION  RATE  (e) 


A  current  modeling  method  for  turbulence  determines  the  turbulent  shear  stress 
from  the  eddy  viscosity  which,  in  turn,  is  obtained  from  the  turbulent  kinetic 
energy  (fe)  and  the  turbulent  dissipation  rate  (e)  from 


\> 


t 


(B.  1) 


where  c^  is  a  constant. 

Values  of  fe  and  e  are  obtained  from  solutions  of  convection  equations  which  are 
partial  differential  equations.  However,  close  to  the  wall  it  has  been  found  that 
values  of  fe  and  £  obtained  from  mixing  lengths  may  be  used  as  inputs  to  the  partial 
differential  equations  for  fe  and  e  with  an  improved  accuracy^  in  the  solutions  of 
the  equations. 

By  equating  the  prodr ;tion  and  dissipation  terms  of  the  fe-equation,  Arora 
et  al.^  related  fe  and  e  to  the  mixing  length  close  to  the  wall  as  follows. 


and 


(B.  2) 


(B.  3) 


Elimination  of  the  velocity  gradient  du/dy  produces  a  direct  relation  between  fe  and 
£  with  %  as  follows.  First  the  relation  for  fe  and  e  are  nondimensionalized  to 


fe  - 


(B.4) 


rr-*” 


num 
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and 


Substitution  of  du  /dy  from  Equation  (10)  produces 


(B.5) 


5 

2 

u 


fe 


21 


*  — i 


l+/l+(2£*)2 


(B.  6) 


and 


v 

4 

u 

T 


e 


i*2 


-3 


l+/l+(21*)2 


(B.7) 


Some  general  properties  of  fe  and  e  are  now  deduced.  At  higher  values  of  1  , 

(/c~/u  )fc  approaches  1  and  (v/u  )e  approaches  1/ H  . 

M  T  4  T  * 

Also  (v/u^.)e  reaches  a  maximum  value  of  1/4  at  £  =  *2. 

This  maximum  value  is  independent  of  the  mixing-length  model. 

* 

In  Figures  B.l  and  B.2,  the  variation  of  fe  and  £  with  y  are  shown  for  selected 
roughness  conditions. 
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Figure  B.l  -  Variation  of  Turbulent  Kinetic  Energy  Near  Wall 


30 


Figure 


APPENDIX  C 

HYPERBOLIC- TANGENT  MODIFICATION  FUNCTION 


A  hyperbolic-tangent  function  may  be  used  as  a  modification  function  for 
arbitrarily  rough  surfaces. 

*  *  2 

M  ■  tanh  (y  /A^)  (C.l) 

* 

The  factor  A.  may  be  correlated  with  the  roughness  function  B  in  accordance 

“  if  it  it  ^ 

with  Equation  (30).  For  A^  *  0,  M  *  1  and  H  -  <y  which  marks  the  limit  of  the 

intermediate  roughness  regime  just  as  for  the  van  Driest  modifications  function. 

*  * 

A  comparison  of  the  integrand  y  (dl^/dy  )  between  the  hyperbolic-tangent 
function  and  the  van  Driest  function  is  shown  in  Figure  C.l. 
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